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Abstract. The singularities of the electromagnetic field are derived to include all the 
point-like multipolcs representing an electric charge and current distribution. Firstly 
derived in the static case, the result is generalized to the dynamic one. We establish a 
simple procedure for passing from the first, to the second case. 



1. Introduction 

In the cases of electrostatic and magnetostatic fields of point-like dipoles, one has 
the well-known procedure of introducing Dirac 5-function terms for obtaining correct 
expressions of the electric and magnetic fields defined on the entire space. The 
corresponding field expressions take the following form [I]: 

Ep(r) = __L p , (r) + = -ipilr) + (E U (!) 

where v = r/r, and 

B m (r) = J^m5{r) + ^ -2 = -fm8{r) + {B) r ^ (2) 

In these equations, by (. . .) r ^o we understand an expression in which the derivatives 
are calculated supposing r ^ 0, representing some well-known expressions of the fields. 
The expressions from equations ([1]) and (j2J) are introduced in Ref. [1] as conditions of 
compatibility with the average value of the electric or magnetic field over a spherical 
domain containing all the charges or currents inside. Another procedure for introducing 
equations (pQ) and ([2D is based on an extension of the derivative didj/(l/r) to the entire 
space [2]: 

a , a i=-f M (,) + i^_%. (3) 

A more pedagogical and suitable approach for understanding the origin of the difference 
between the electric and magnetic cases, is done in Ref. [3]. Refs. jl] and [5] contain 
generalizations of the equations (Jl]) and (J2J) to the dynamic case for oscillating electric 
and magnetic dipoles. 
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The objective of the present paper is to establish the singularities of the 
electromagnetic field associated to a system of electric charges and currents assimilated 
with a point-like multipolar system. These singularities are established for an arbitrary 
multipolar order for both, the static and the dynamic cases. 

In section [2j the procedure of separating the 5-type singularities of the wave 
equation retarded solution is described. In the next section, based on the exterior 
solution of a given static electric charge and current distribution expressed in terms 
of electric and magnetic multipoles of this distribution, i.e. the multipole expansion of 
the electromagnetic field, we separate the 5- type singularities of the fields E and B. 
The calculation is based on the invariance of the static electromagnetic field to the 
substitution of the Cartesian tensors representing the multipole electric and magnetic 
moments by the corresponding symmetric trace free (STF) projections. In section HJ 
the singularities in the dynamic case are established by employing also the invariance of 
the electromagnetic field to the substitution of multipole Cartesian moments by STF 
tensors, but this time, generally different from the corresponding STF projections. 

In Appendix A, we give some furmulae for the lower order spatial derivatives of the 
spherical wave retarded solution used in the current exposition. In the second appendix, 
the reduction to the STF tensors in the dynamic case and for low orders is presented. 

We point out that the formalism presented in this issue has as mathematical basis 
the properties of the irreducible tensorial representations of the proper rotations group 
[6]. Although, for pedagogical and larger accessibility reasons, we give an explicit 
calculation based on the properties of the tensor contractions such that the procedure 
has a simple algebraic character. A counterpart of the procedure used in the present 
paper could be represented by the technique of the spherical function expansions and 
the cited issues can be a basis for such an approach. 

2. Some delta-function identities 

The treatment of some delta-function identities in Ref. [2] can be easily generalized to 
obtain the necessary delta-function identities in the dynamic case. Let us consider the 
solution f{t — r/c)/r of the wave equation. 



In the points different from the origin O, this function is a solution of the homogeneous 
wave equation. The multiple partial derivatives of this function are given, for r ^ 0, by 
the formula: 



where r = t — r/c. The coefficients are symmetric in z'i, . . . , i n and can be expressed as 




(4) 





{h h 



...5; 



(6) 



k=0 
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In the last equation, [/?] is the integer part of (3 and v = r/r. By i; ilcJn i we 
understand the sum over all the permutations of the symbols i q giving distinct terms. 
For the objective of the present paper, the coefficients D in equation ([6]) can be 
calculated directly from the successive derivative operations. We give in appendix A 
the corresponding expressions for n < 4. 

From equation (j3J), it is obvious that we have to consider the derivatives of f(r)/r 
extended over the entire space as distributions (generalized functions) which contain 
singular distributions, as ^-functions for example, as separate terms. Let be a function 
F(r, t) and suppose the existence of the integral of the product F(r, t) 0(r), with 0(r) 
an arbitrary smooth function (a test function from the domain of the distributions), on 
the spherical region T>r, with arbitrary radius R, delimited by the spherical surface E# 
with the center in O. 



d 3 xF(r,t) <p(r) = lim 



d 3 xF(r,£)0(r) 



(7) 



J £< r<R 

The integral can be expressed excluding from the domain T>r a spherical domain of 
radius e centered in O. Writing this last limit of integrals, we can interpret the function 
F(r, t) as a distribution defined by 



(F(r,t)) r ^ , 0(r) 



lim 



V R \V S 



d 3 x (F(r,t)) r ^<f>{r) 



This distribution can be extended such that its support includes the point O. A new 
term 9(e — r) F(r, t) can be naturally introduced by the identity 

F e (r, t) = die - r) F(r, t) + 9{r - e)F(r, t), 

associated with the extension of the integral to the entire domain T>r\ 



{F(r,t),<f>(r)) 



lim 



/ d 3 xF(r,t)0(r) + / d 3 xF(r,t) 

JV e JT> R \D e 



Moreover, we suppose the existence of the integral (171) for the partial derivatives of 
F. Let us consider the partial derivative diF(r, t) and the problem of extending this 
function as a distribution. The definition (181) becomes: 



{F(r,t),<j>(r)) 



lim 



/ dSviF(r,tU(r) - f d 3 xF(r,t)d, 
f d 3 xdiF(r,t)(j)(r) 



(9) 



where S £ is the sphere of radius e centered in O and the Gauss theorem was employed. 
Let us apply this definition to the derivative didj(f(r)/r) and, for simplifying the 
notation, let 

f(r) 



D h...i n (r,t) 
such that we can write 



d h . . . ^ 



(D 



lim 

£->0 



d5* V; d. 



f(r) 



d 6 xd. 



fir) 



0, 



d 3 x (aa^ 

v R \v E r 



(10) 
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Since the last integral on the domain V R \ V £ represents the distributions associated 
with the F- expressions for r ^ and, in the case of the electromagnetic field, will be the 
well-known expressions of the multipole expansions, in the following we consider only 
that part of (A/) containing singular distributions with point-like support i.e., actually, 
the difference 

f{r) 

e-*oJv R \v s v " J r 



(Aj) (0) , 0) = (A,, 4>) - Km /_ d\(^^)0(r) 



lim 



dS Vi dJ-^-cb(r) - ( d 3 xdi^dMr) 



s £ r Jv e r 



ill) 



By Dt TQ \ we denote a distribution having as support the point given by the vector r . 
In the first limit, we make use of equation (|3j) and of the expression of the coefficients 
Cj- from Appendix A: 

l im J dSV i d i ^0(r) = -lim/ dS v iVj [-f{j) + \f(r)} cf>{r) . 

Inserting the Taylor series of the function 0(r) and since on the sphere r = e, we can 
write 



lim/ dSuidi^ 



= -4vr lim J dQ(u) v iVj f(r e ) + /(r e )] [0(0) + ev k (<9 fc 0) o + . . .} (12) 
where r £ = t — e / ' c. Let us introduce the angular average: 

<^)> = ^/ g(")M{v). (13) 

Particularly, we have the well-known formula [TJ: 

(u v-)-[ °' n = 2k+l, 

(n+l)\\ %1*2 • • • ^i n -l»n}> H = ^ fc = 0, 1, . . . 

Excepting the term containing the product /(r e ) 0(0), all the terms in equation (1121) are 
proportional to positive powers of e and, consequently, vanish with e — > such that 

lim I dS Vi 6>j — 0(r) = - 4tt ( z^) /(t) 0(0) 
Considering the second integral in the right-hand side of equation ( fTTT) . we can write 



d 3 x d,— 9i0(r) 



= - lim [r 2 dr [ dn(u) u,[-f(r) + -f(r)][(d l( p) + ru^dA^o + . . .] =0 (15) 
such that finally 

($3,-—) =~Wt)8(r). (16) 
V r /(o) 6 
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In the static case, f(t) = 1, equation (ITSl) becomes equation (3) from Ref. [2] after 
adding the regular distribution represented by the last integral from equation ffTUj) . 

Let us consider the distribution Dij k . Considering only the part having O as 
support, 



((Aifc) (0) , 0) = Jim £ dSuidjdk^p-^r)- d?x {3,3^-^)8, 
The surface integral becomes: 



.(17) 



/ /( r ) / 1 •■ 1 

lim & dS Ui djd k <f>(r) = lim * dS Vi\-^-VjV k f(r) H -(3^^ fc - 5 jk )f(r) 

£^oJs e r e^oJs e err 

1 



+ ^( 3 ^ - 5jk)f(r)]<l>(r) 
and, introducing the Taylor series for 0(r), 

lim / dS i/j djd k ^-^- <p(r) = Air lim / ^-ViVjV k f(r £ ) 

+ -(3^^-^ - Vi5 jk j(T E ) + -(3^-z^ - Vi5jk)f{r E ) [0(0) + ez/j (<^0) o + . . .]\ . (18) 

All the terms containing the products /(r e )0(O) and /(r e )0(O) give null results being 
multiplied by averages of three factors v. The other terms containing the factors /(r e ) 
and /(t e ) give also null results being proportional to positive powers of e. The term 
provided by the product /(r e )0(O) contains the factor but it has also a null limit 
since (3viVjVk — ViSjt) = 0. Only the product /(r £ ) (<9z0) o has a factor independent of e 
having a limit different from zero. 



Therefore, 



lim / dSvsdidk^ 



1 „ „ 1 



4vr (-%4o - 3^5 jfc ) /(*) (<9,0) o (19) 



Concerning the integral on T> £ from equation ( fPTl) . we have to observe that, 
beginning from this derivative order, there is a non-zero contribution for e — > 
[2]. Indeed, introducing equation (fT6l) in equation (fT7j) and observing that the term 
(djdk(f(r)/r)) r ^ Q gives a null contribution to the limit for e —> 0, we can write 

-lim / d\d 3 dj-^ d i( f>(r) = %\ d 3 x/(t)^(r)^0(r) = ^/(t)^0) o . (20) 
e^oJv E r 3 Jv E 3 



Finally, equations (ITTj) . (1181) and (jl9j) give 

47T 

(A jfc ) (0 ) = -y /(*) S w d k} S(r) . (21) 

For the static case, f(t) = 1; this result corresponds to equation (4) from Ref. [2]. 

Obviously, this procedure becomes very complicated for higher order derivatives. 
Fortunately, for the electromagnetic field some invariance properties allow a considerable 
simplification of such calculations. 
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3. Singularities of the electromagnetic field: the static case 

Let us consider the multipole expansions of the electrostatic and magnetostatic fields. 
Given an electric charge and current distribution with support included in the domain 
T>, the scalar potential is expressed in the exterior of this domain by the following 
multipolar series: 

In this expansion, the coordinate system origin O is supposed in T> and is the n-th 
order electric multipolar moment defined by the Cartesian components in the general 
dynamic case: 

P h ... in (t) = J v d 3 xx H ...x ln p(r,t): P< n \t) = J d 3 x r n p(r, t). (23) 

v 

In equation (1221) we employed the following notation for tensorial contractions: 

{Ail—in-mjl—jm-Bji—jm i n > m 
A h .., n B h .., n ,n = m . (24) 

Aj 1 ---j n Bj 1 ...j n i 1 ...i m _ n , n < m 
For the vector potential in the exterior of the domain V, we have 

71-1 



A (r) = — E , V x (V^HMW 



= E • • • d^^f^, (25) 

where M*"^ is the magnetic n-th order moment defined by the Cartesian components 
0: 

M h...i n (t) = ——r / d 3 « 11 ...i i „_ 1 (rxJ(r,t)). ) (26) 
n + 1 Jp 

or, with tensorial notation: 

M (n) (t) = — ^— / d 3 x r n x J(r,t) . 
v 

For the multipole expansion of the electric field E(r) = — V$(r), equation (1221) 
becomes: 

1 ( — -])n-l p(n) 1 f_1\«-l p. . 

^(r) = T~Y. V n+1 || — = -i-e, £ { -^d i d il . .AKP^.. (27) 

47r£ o^l n\ r A-keq ^ n! r 

where, for simplicity, the electric charged system is considered neutral (Q = 0). 

The corresponding expansion of the magnetic field B(r) = V x A(r) is given by 

^)->x E Hp V xfv-||^ 
4vr & n! V " *" / 



HQ ^ (-1 



in-1 



E 



47r n! 



V • V n || - A V 



(28) 
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We have to search the singularities of E and B given by equations (1271) and (12"8"|) . It 
appears that cumbersome calculations are involved for higher n if we apply the formulae 
for derivatives of an arbitrary function /(r)/r as in the previous section. However, we 
can employ an invariance property of the electromagnetic static field to the substitutions 
of all moments P^ n ^ and for all n by their corresponding symmetric and trace- free 
STF projections "PW and Ai^ [9], [TDl E]. Retaining the notation p and m for the 
first order moments, this invariance stands for the invariance of the electromagnetic 
field, in the static case, to the following substitutions: 

p, P( 2 ), p( 3 >, ... - p, V {2 \V^... ; 

m, M (2) , M (3) , ... -> m, Al (2) , M (3) .... (29) 
These STF tensors can be expressed by the following formulae: 



(2n-l)!!7o rv 7 ^•••^ r 

/ i \n+l n « i 

M u i =7 -rr ttt F / d 3 xr 2n+1 (J x V), ~, (30) 

1 " (n + l){2n - 1)11 ^Jv K hx ll ~ l " r ' v ' 

which, actually, differ from the projections by numerical factors. Here, 

8 W = 8 8- 8- 8 
u n...i n u n ■ ■ ■ u i\-i u i\+t • • • u i n ■ 

Let us consider the delta-singularity corresponding to the electric dipolar field: 

a result obtained applying directly equation (TIB]) for f(t) = pj. Concerning the 4-polar 
term from E, we consider firstly the expansion (|27[) expressed by the primitive moments 

p(n). 

E^(r) = --^-e l 8 l d J d k — (32) 

87T£:o r 

and from equation ff2TT) : 

= ~^T e i PjkS{ijd k} 5(r) = y^r^ P jk (5 i: jd k 5(r) + 5i k dj5(r) + 5 jk 8i5(r)) 

= J-ei (2P ij d j 6(r) + P^-^^r)) = ^-e t (Pyfi^r) + ^^(r)) . (33) 

Note that P^ 2 ^ is symmetric. 

The substitution Py — »■ V%j in equation (1331) gives 

JSSS = r^7^(r) = -^p( 2 )| | V*(r), (34) 



since T^ 2 " 1 is trace-less ('Pjj = 0). 

From this simple example, it becomes obvious that it is more simple providing 
formulae for the delta-singularities directly for STF tensors. Let us consider the 4-pole 
approximation using the definition 



w I 8ire e^o Jt> £ \ r 



-e, lim 



87T£ £ -*° 



I dSV jk Ui (djd k -} <p{r) - f d 3 xV jk djd k - d? 



(35) 
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In the last integral over T> £ from equation (1351) . we can introduce equation (1161) obtaining 



lim 



<TxV jk 



Air 



—5 jk 5(r) + [djdk- 



r J r^O 



di<f)(r) = lim f d 3 xV jk (djd k ^) <%0(r 



lim / d 3 x V jk \ (3ujU k - 5 jk ) d^r) = 3 lim / d 3 x \v jk VjV k <9;0(r) 



3 lim 



f dVtV jk VjV k [(di(f))o + rui{didi4>)o + . . .] 



0. 



since the term Vj k VjV k (di(f))Q provides the expression Vj k 5j k = and for the other terms 
from the series has null limits containing, after integration over r, positive powers 
of e. The same result can be obtained by a repetition of the integration by parts in the 
last volume integral. Therefore, we can write 



eilim/ dSV jk v,i (djd k -) <p(r) 



1 1 

- — ei \imV jk - {vii^VjV k - 6jk) (0(0) + eui(di4>) + . . .))) 



1 



- TT-ei'Pjk (ZViVjVkVl) (di<j))) 
ZEq 

_ J_ P (2)||( V( 

5e 



1 



I0e 



V jk {5{ij5ki}) {.di(f))o 



i.e. the result 

Let us consider the general terms from the series representing the electric and 
magnetic fields in equations ( 1271) and (1251) with the substitutions ( 1291) performed. 
Denoting by T^ n ^ an n-th order STF tensor, we observe the presence in both equations 
of an expression of the form 



Bidi d h ... d in 



T; 



Searching the extension of this function to the entire space, we introduce the distribution 
I(TW) by 

T(")\ 



J(T^), 



Writing the last equation as 



lim 

e^O Jj), 



(n)^ 



lim 

e->0 



dSl/ V n | 



d*x V V 7 



d 3 x V n | 



V</>(r] 



(36) 



(37) 



we can consider separately the limit of the surface integral: 



iim£ds,(vii^) m 

Bi lim / d5 vJit.-.i-n (dh •••<%„- 



(r). 



Applying equation (JHJ) and the Taylor series of </>(r 



47re,- lim ; — T,-. ,■ i.Vi C\ n,n } Vi , , 

o=0 



X I • • • din+a 



(3? 
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Let us evaluate the tensorial contraction which is present in the general term from the 
previous series: 

(39) 

From equation ([6]), we can easily see that all the terms containing at least a symbol 
5i q i a with 1 < q, s < n give null results by contraction with T^ n \ An exception is 
represented by the term containing the factor Dq 1 '^ v% x . . . Vi n which, together with the 
factor za ,-. . . . Vi , ia-, can give results different from zero. Since 

D^ n) = (-l) n (2n-l)H , (40) 

the series ( |38l) can be written as 

oo a—n+1 

(l a (T^),<f>) =47rc i (-l) n (2n-l)!! Urn £ — 

x Ji,...^ (v h ... v in u in+1 . . . v in+a (d in+1 . . . d in+a 0) Q . (41) 

Firstly, we have to see what happens for negative powers of e in the terms of the series 
from the last equation. These negative powers are factors in the terms verifying the 
inequality a < n — 1. The expression is different from zero if and only if it contains 
some terms in which all the indices i\ . . . i n are present as index pairs of Kronecker 
symbols 6^ or 5a q with j, q — 1, . . . , n and k £ {n + 1, . . . , n + a}. Otherwise, one 
obtains null values because of the contraction with an STF tensor. For a < n — 1 such 
terms cannot exist. 

For a > n — 1, the corresponding terms from the series in equation (jUJ) contain 
positive powers of e. Consequently, the corresponding limits for e — > vanish. From 
this series only the term for which 

a = n - 1 (42) 

can be different from zero and the result of the limit in equation (1411) is given by 

47r(-l) n (2n - 1)!! 



Irr(T [nY 



(n- 1)! 

x e.i T h _.. in ( u h . . . v in u in+1 . . . v^ n _ x (<9 Wl ...di 







47r(-ir(2n-l)!! / (n) 2n 
(n-l)\ \ 11 



jn—l. 



(43) 



TTl— 1 



5(r) . (44) 



The corresponding 5-function is therefore 

(n 

Let us consider the contraction 

( T W|| Z , 2 "') = T n ... in { v ix ... v in v h ...v in ) . 

To this contraction contribute only the terms from the average of the //-product not 

containing factors S ikil , with 1 < k, I < n. According to equation f|T4|) . the terms giving 

non-zero contributions are of the form 
1 



(2n 



6i 1 j 1 . . . Oi n j n 
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and all such terms are obtained considering the n\ permutations of the indices j\. . . j n 
in this product. Therefore, the final expression in equation fj44H is given by 



I*(J (n) ) = -^l^ n) W^ n ~ X Kr) ■ (45) 
Let us evaluate now the volume integral in equation ([3 



J(T (n) ), V<j)) = \imJ^ d 3 x fv n | 



Vcb(r) 



(46) 



r ( T^ n )\ r ( T( n )\ 

= lim i dSv\\ V" _1 || V0(r)-/ d 3 i V n ^|| V 2 ^ 

e^O [Js e \ r J JVe \ r J 

The part corresponding to the surface integral can be written as 

■UT {n) ), V0) = ej tl ... tn lim£ dSX (\ . . . c^-) <9 4 0(r) . (47) 
Introducing the Taylor series for 4>(r) and standing out the average over is, we obtain 

oo a— n -\-2 

J CT (TW), V0) =4vreaim E— ^-T^ (cfc^V ■ ■ ^w- 



x ( <% n+a • • • di n+a d4 ) n • 



Analogously to the reasoning from the previous case, we can see that the limit is zero 
since for a = n — 2 



Y (re— l,n— 1) 

1 i\...i n \ <~ 

Therefore 



1 ti...i„ \ ^i x ...i n -\ V in V in+\ ■ ■ ■ v ii n -ij ~ U - 



J(TW), V0) = -lim / d 3 x V"" 1 !! V 2 0(r) , (49) 



and by repeatedly applying the procedure, all terms cancel. Finally, only the surface 
integral from equation (157)) gives a limit different from zero and 

I(T (n) ) = 7 ff (T (n) ) . 

Introducing the result (1431) in equation (127)) . we get 



^- (n-DuLl)^ "' 11 ^"^- (50) 



With this formula one obtains easily the results for n — 1 and n = 2. If we take the 
3-rd order term of the electric field, equation ( 150)) gets the form: 



= -^P (3) ||V 2 5(r) = - J-c^WCr) . (51) 

14 fcQ 14 fcQ 



For the 4-th order multipolar term, equation (150)) reads 

,(4) _ _J 

'(°)"54e 

and so on. 



^j = ^p (4) ||V 3 5(r) (52) 
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Concerning the magnetic field, selecting the singular point-like part of B^ from 
equation fl28l) . we can write 



B 



(n) 
(0)' 



^4^hm I 



7TI+1 1 



rn-1 



A- 



:-i) n -Vo 

47T77,! 



lim /d 3 xM (n) || V n+1 - + 47r(v n ~ 1 5(r 

e^Oj r \ 



4:71 n\ 



lim / d 3 x ( M (n) ||V n+1 - ) 0(r) + ^.M (n) 
e^Oj V rJ TV. 



jn—l. 



since A(l/r) 



Air8{r). Written explicitly, this last equation becomes: 
B&l <A = ( ~ 1)W l^ e t lim / d\rM r ...,j), (d n . . . d in -) J>(r) 



'(0)' 



Mo 



+ -jeiA^n... in _ 1 i(9 il ...9 in _ 1 



(53) 



The integral on T> £ from equation (1531) is of the same type as the integral from equation 
(|36|) and we can apply the result (145]) obtained in the electric case, such that 

' (-l)>o (-1^" 



B 



(0) 



i.e. 



B 



(») 



(n-l)!(2n + l) 
-lf-Vo (n+1) 



+ 



Mo 



?>,! 



M (n) ||V n - 1 <J(r). 



(0) " n! (2n + 1) 

Let us apply equation ( 1541) for lower multipolar orders. 
For the dipolar case: 



(54) 



For n = 2, the quadrupolar case, 



2/ic 



(2) 
(0) 



For n 



B 



(3) 
(0) 



*!*MV>\\V8{t) . 
10 11 K ' 



(55) 



(56) 



(57) 



Before passing to the problem of the variable electromagnetic field, we point out 
some simple consequences of the formalism employed in the present section concerning 
the integrals over a spherical domain T>r of the electrostatic and magnetostatic fields, 
supposing the support of electric charges or currents included in Vr. 

Using equation (j22|) for representing the electric potential on the spherical surface 
we can write, employing the invariance of E to the substitutions (1291) . [1] 



drxE(r) 



dSv$(r) = 



E 



\n-l 



dS v V" 
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By applying the Gauss theorem in the above integral, there were supposed fulfilled all 
the conditions of continuity of the field on the surface S of the domain T> where the 
electric charge is distributed. 

Considering the surface integral corresponding to the n-th order, 



■R 



dS v V T 



An 



R 



n-l 



Bi ( Vi C) 



(n,n) 



i ^ti...t r 



V 



Excepting the case n = 1, the contractions of the angular averages with the STF tensor 
?W give null result such that 



and 



An 

I R = Anei( ViVj) pj = — —p 



f d 3 xE(r) = -^-p 



i.e. equation (4.18) from Ref. [T]. 

In the case of the magnetic field, using the exterior solution for A given by equation 
25]), we can write 

fi ^ (-IT- 1 



d^xB r) 



dSu x A 



E 

n>l 



111 



dSux 



V x V 



r 



Considering the surface integral corresponding to the n-th order, 
J R = eiE i:jk e k i q dS h>j (did h . . . di n ~) Mi u .. in 
dQ(u) 



i£ijk£klq V 3 °«...i„-i l^il-i-n.-! q 



An 



As in the electrostatic case, excepting the case n = 1, the above contraction gives null 
result such that 

J R = -Ane ijk e k i q ( VjV\) m q = — ^ (SuS jq - 5 iq 5ji) 5jim q = 



and consequently, 



d 4 xB(r) 



2fi 



iv R 3 
i.e. equation (5.62) from Ref. [I]. 



rn 



4. Singularities of the electromagnetic field: the dynamic case 

In the dynamic case, the multipole expansions of the potentials in an arbitrary point 
exterior to the domain D are given by the expressions [TT] : 

4vre n >! n\ r 
1)n -i 



n>l 



111 



V x (V 



n-iM^(r) 



(58) 



Singular behaviour of the electromagnetic field 



13 



As in the static case, for simplicity, we consider a neutral electric system (Q = 0). We 
write the multipole expansion of the electric field E(r, t) derived from the potentials 
given in equation ( 1581) : 



E(r, t) 



Aire 



n>l 



VfV 



n„ P(n) W- 



IV. 



n-l 



_Ivx[V-||^ 



The multipolar expansion of the magnetic field is given by 

.l)n-l 



B(r,t) 
Since 



47T 



n>l 



»! 



v x (V"- 1 !! - 1 J ) + v 



n-l| 



V x 



V x V 



n-l| 



M( ft )( 



MW(i 



A V 



n-l 



(59) 



(60) 



T 



equation (1601 can be written as 



B(r,t) 



47T 



E 

>i 



-i 



,n-l 



'/?! 



M( ft )( 



A V 



n-l 



n>l 



-1 



,n-l 



-V x V 



n-l 



P T) 



(61) 



Equations ([BUI) and flBTT) represent the fields £J and S for r 7^ 0, where the multipole 
tensors considered at t — r and divided by r are solutions of the homogeneous wave 
equation. This property can be considered in any processing of the field multipole 
expansions. The corresponding delta-singularities, i.e. including the point O in the 
domain of the field, will be expressed searching the extensions as generalized functions 
of the final expressions. 

As in the static case, for higher orders of multipolar terms, n > 2, the labor involved 
in the calculation becomes prohibitive if we operate with the primitive moments 
and M( n ). Fortunately, there is an invariance property of the electromagnetic field in 
the dynamic case, too, which allows us to replace in the expressions of E and B the 
tensors P < - n - ) and M*™) by symmetric and trace-free (STF) tensors for any n: 



p(n) 



M (n) 



M 



CO 



(62) 



~( n ) ~(n) 

P and M are STF tensors (see Refs. pH [131 E] for n < 3 and [LTJE5] for arbitrary 
n). Theoretically, each such STF tensor is an infinite series, but practically, we have 
to calculate only a finite number of terms from the multipole expansion such that the 
respective tensors are represented by finite sums. Consequently, we can operate in the 

^ ~ (2) ~ (3) — -(2) 

following with the STF tensors : p, P , P , . . . for the electric cases, and m, M , . . . 
for the magnetic ones. However, for not invoking directly some results from the issues 
cited above, we introduce in Appendix B|the lower orders STF tensors employed in the 
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present paper, together with the introduction of the STF projections of the primitive 
momenta. In every case, we will establish what are the changes of different momentum 
tensors which compensate the effects of such substitutions. 

Let us consider the problem of expressing the singularities of the electromagnetic 
field in O when we are interested only in the contribution of the electric dipolar moment 
p(t). In this case we have an elementary charged system characterized only by the 
dipolar electric moment as, for example, the harmonic oscillator. Another case is when 
we are interested only in the first multipole approximation for a complex system, as an 
atom. For n = 1, considering the contributions of the electric dipolar moment, equations 
((59]) and flBJ give 



-i7-p(t)S(r) . 

6Eq 



(63) 



The singularity in the case of a magnetic dipole is 



r»(l,m) 
^(0) 



i.e. 



An 



didj 



mj {T) 



771,; T 



-mj(t)SijS(r) +mi(t)S(r) 



B 



(l,m) 



2^0 



m{t)5(r) 



(64) 

from the static case, one 



'(0) " 3 

Comparing equations (1631) and fl64"|) with equations ([I]) and 
observes that the results for the dynamic case are obtained from the static one by the 
substitutions p —> p(t), m — > m(t). 

In equations ( 1591) and (1611 . for the multipole expansions of the fields E and B we 
have a first general term of the form 



7TI+1 



T 



eA d n ... d in 



T; 



(r) 



introduced for r ^ and with a STF tensor . Let us consider the extension of 
this function as a distribution isolating the point-like singularities by the definition 

,T (n) (r)\ ,, , 



l(T (n) (f) 



lim 



ej lim 



d 6 x V 



n+l I 



Jn 



d x [didi,... d in 



(65) 



The partial integration becomes 



!(TW(f)), 



lim 



dS v V r 



d^x V T 



||T (n) (r) 

r 



(66) 



Considering only the surface integral from the last equation, employing equation ([5]) 
and the Taylor series of the function 0(r), after a regrouping of factors, we obtain 

n oo c a— l+l 



X UiC, 



1=0 a=0 
(n,l) d n 'T^...^ 



(r«) 



dt 



n—l 



Pi 



(67) 
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where T e = r — e/c. As in the static case, we can write this last expression in terms of 
the averages upon the directions: 

n oo a— l+l An—l 

(I„ (TW(*)) , 0) = toetVmE £ (cfclT n ..,^K +1 . . .^ +Q „<) 

1=0q=0 L Ut 

x (^ +1 ...9 Jn+Q 0) o . (68) 

The general form of the coefficients C^ n '^ is given by equation ([6]) and one can notice 
that only the term D^' 1 ^ z/^ . . . gives a non-zero contribution to the contraction with 
the STF tensor T^ n ^. Therefore, we can write 

n oo a— l+l 

I. (T<»>(i)) , 4.) = 4,e,lim i: £ 

X^Z7 T U...^( r £ ) ( fil • • • ^in • • • ^n+a "i) ■ • ■ 9 in+a </>) Q (69) 

In this equation, the factor e a ~ l+1 represents a negative power of e for a < I — 1. 
Since Z < n, we can write 

a < n - 1 . (70) 

In this case 

(u h . . . v in u in+1 . . . u in+a ViJ ^ =7 T il ... in (r e ) = 0, (71) 

since n + a + 1 < 2n. Finally, the combinations of the factors e a ~ l+1 for a < n — 1 
with the terms from the Taylor series of T^.. i n (r £ ), 

00 (-1) A £: A d A 

can give non negative powers of e, but equation (1711) remains valid. 
The terms not containing e as factor verify the equality 

a = I - 1 (72) 

and they can give contribution different from zero. 

For a > I — 1 the corresponding terms contain positive powers of e and have null 
limits for e — > 0. Consequently, the sum over a is limited at a = I — 1. 

Summarising, one can write the following expression for the limit in equation (166]) : 

/jT<">(t)),*)= 4™, 



1=1 a=0 L " 



d n-l 



) (^n+l---5 in+Q 0) o .(73) 



dt n - 

For Z < n, equation (1721) becomes a < n — 1 and equation (17T1) is verified; in the sum 
from equation (1751) remain only the terms corresponding to I = n and a = n — 1 which 
give: 



x T^.^Xt) . . . ^ in+1 . . . i/<) (<9 in+1 . . . d i2n _ x 4>) q . (74) 
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Comparing this result with the corresponding distribution (1381) from the static case, 
one can see that equation (174"1) can be obtained from equation (I43p by the substitution 
j(n) _^ TM(t). 

The limit of the volume integral from equation (|65|) can be written as 
J (T {n) ) , V</>) = lim d 3 x ^V n T(n ^ T ^ V0(r) (75) 



lim <^ & dS 



v V r 



|,T (n) (r)y 




V0(r) - J^d 3 x 


r J 





|V 2 0(r) . 



The limit of the surface integral from the last equation can be processed as in the case 
of equation (l67j) . by writing 



ei dS ( dh ■ ■ ■ d ^-i^^ 



d 

x / c!-™~l ,n ,Vi Vn . . . Vi , \ r^-Tj- ; (t>) f Aft . . . di , 0) 



/• 



^)^(r) = eaim4/EE L ^T 



n—1 00 ( 1 ^A g-"— '+1 



df 



(76) 



Only the terms with a — / + 1 = can give contributions different from zero but, since 
I < n — 1, a verifies the inequality a < n — 2. Consequently, the average from the 
last equation contains at most a symbol with j, k < n and the contraction of this 
average with the STF tensor is zero. Finally, by reccurence, one can verify that 
the limit ( 1751) is zero such that equation ( 1741) is the result for the limit from equation 
(J65])- 

Similar considerations lead to the conclusion that the second term, proportional to 
V n_1 ||(P (r)/r), in the right parenthesis from the expression (|59l) of the electric field 
gives null contribution to the point-like singularities. The contribution of the magnetic 
moments to the point-like singularities of the electric field are defined by a distribution 
of the type 

Tji...i n _i k( T ) \ ,/\ 



ei lim e ijk / d x \dj d h . . . ■ 



e, lim 



s dS uj ld h . ..di n ^ 



d 3 x \d il . . . 



T 



% 1 ...i r 



-At)' 



d j( f)(r) 



(77) 



iv E \ r 

Concerning the surface integral from equation ( ITTj) . introducing, as in previous cases, 
equation ([5]) for partial derivatives and the Taylor series for 4>(r) (with the same notation 
for the various parameters), we notice that each term contains the contraction 

.a-l+l _ / £j(n-l,n-l) 



This contraction can give contributions different from zero only for a = I — 1 < n — 2. 
As one can easily see, for such values of a this contraction vanishes either because of 
the presence of the symbols Si n i with q < n — 1 or of 5i ; with q, s < n—1. The same 
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conclusion applies to the contributions of the electric moments to the magnetic field as 
seen from equation (16T1) . 

The last type of singular distribution we have to search is for the term 
A[V n ~ 1 ||(T (n )('r)/r)] from the expansion (}6Tj) of the magnetic field. The corresponding 
point-like singular distribution can be obtained from the definition 

(Q(TW), <f) =e k fmj^xd il ...d in _ 1 Llil^lM\ 0(r) 

= -47rlim/ d 3 x J h ...i n ^ kit) <p(r) d h . . . d in _ x S(r) 

= -4 7 rT il .. in _ lfe (t) ( d h . . . d in _, 5, <j>) • (78) 

Equations (1711) and (1751) . inserted in equation (loTi) . lead to the same conclusion as in 
the electric case: the distribution B^} in the dynamic case can be obtained from the 
static case, given by equation (1541) . with the substitution — > Ad^ n \t). 



5. Conclusion 



Based on the results of the previous section, it is obvious that in the expansions 
and (16T1) of E(r,t) and B(r,t), after performing the substitutions fl62l) . only the terms 




r 

for the electric field, and 



V|V n || ^ | , A|V n 

r 

for the magnetic one give contributions to the ^-singularities of the electromagnetic field. 
These contributions can be obtained by simply performing in the corresponding static 
expressions ( l50i) and ( |54|) of E and B, the substitutions 

V (n) ^ p^( t ) ; M (n) ^ n \ t) (79) 



In Appendix B the reduced STF moments for the first multipoles are established. 
Considering the substitutions (179!) performed for all the multipoles, obviously P (t) 

— (n) 

and M (t) are actually represented by infinite series. Practically, as seen from this 
appendix, for a given multipole approximation, one deals only with finite sums. If, for 
example, we are interested only in the contributions to the singularities of the field from 
the electric octupole (n = 3) and magnetic quadrupole (n — 2), besides the contributions 
of the STF projections T*^ and M.^ 2 \ contributions given by equations (15T|) and (156|) . 
we have to consider also a contribution of the type (I3T1) of the electric dipolar toroidal 
moment T defined by equation (IB. 201) . This last contribution is given by equation (I3T1) 
with p — > p = p — T/c 2 . If we are interested only in the contribution of the magnetic 
quadrupole, besides the contribution of the corresponding STF projection, we have to 
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consider the contribution of an electric dipolar moment A'p given by equation (IB. 16|) . 
and so on. 

We point out that there are alternative ways for introducing the moments P and 
M. Obviously, the results from the present paper obtained through the formalism of 
the Cartesian tensors, can be obtained employing the expansions of the electromagnetic 
field in the formalism of the spherical tensors as in Refs. [1] and especially [12], [13] . One 
can employ the multipole moments derived by this formalism and then introduce the 
connection between spherical and Cartesian components of the multipoles [H] . 

The results of the present paper, if correct, could be useful in studying some 
problems of interest from the atomic, nuclear, elementary particles, and even gravitation 
physics. 

Appendix A. Expressions of the coefficients C^ n >^ 

The coefficients C from equation ([5]) for < n < 4 are: 
= 1 ; 

-v h C\ X) = -v % - (A.l) 

UiUj, Cj 1 = ZitiVj - 5ij, C\f 2) = 3^ - Sij ; (A. 2) 
= -ViVjV k , C§£ = -6h>iVjV k + 8 {ij u k} , 

-- - lBi/ii/jVk + 3 8{ijV k }, C (3,3) = - VaViVjV k + 3 %z/ fc} ; (A. 3) 

ViVjV k V U 

\§ViVjV k vi - 5 {ij v k v l} , 
45uiUjU k ui - 65 {ij v k v l} + 
105 ViVjV h vi - 15 5 {ij v k v l} + 3 5 {ij 5 kl} , 

105 v^^i - 15 8{ijV k vi} + 3 8{ijS M y. (A.4) 



C 



(1,0 



a 



(2,0) 



r»(3,o) 

(7(3,2) 



r»(4,o) 
r»(4,i) 

^ijkl 

r»(4.s) 

^ijkl 
(7(4,4) 



Appendix B. Reduction of the multipole tensors 

Let us consider the electric quadripolar and magnetic dipolar contributions to the 
electromagnetic fields. The 4-polar electric moment P*- 2 -' represents the first tensor 
for which one has to apply the substitution by the corresponding STF one. The tensor 
P^ 2 ) is symmetric and the corresponding trace-free projection, denoted by can be 
established as a combination of the form 

V,j P, 7 A8; r (B.l) 

The parameter A is determined from the condition Vjj = 0: 

A = ^P 4l . (B.2) 
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Firstly, we observe that E and B are not changed by the substitution P( 2 ) — > 7*^. 
Indeed, equation fl55l) reads: 



p(2) _,-I>(2) 

E(r,t) P ^ E(r,t) 



8ne 



A(r) 



I a .A(r) 
c r 



(B-3) 



Since r ^ 0, equation (1B.3|) shows the invariance of the electric field to the considered 

substitution . Considering the change in the expression floT|) by the substitution 
p(2) _> p(2) > we 

A(r) 



p(2)^p(2) , ,\ oor A(t) 

B(r,t) — ► B(r,t) + —eieijkdjdidik—^- 



8n r 



B being also invariant to this substitution. Let us consider the substitution 

p(3) _> p(3) 



0, (B.4) 



(B.5) 



in the multipole expansions of the fields E and B. We have to express the components 
Vijk in terms of the components P^. Obviously, the relations between the components 
of the tensors and T 9 ^ are of the following form: 



Considering the conditions on the traces Vuj we write 



(B.6) 



-p.. 



1 1 J 



5 Aj = 0, 



i.e. 



A, 



1. 



1 



5 Jv 



d 6 xr 2 Xi p . 

The effect of the substitution (IB.5P in the expansion (1531) of E is given by 

A,fr) Ai (r) 



(B.7) 



E^(r,t) -+E®(r,t) 



24vre 
1 



r 



247T£ C 2 



Si did f^ 3 



set- 



or with the tensorial notation, the change of the field E is given by 



E(r,t) —5 E(r,t) 



2Ane c 2 



V V • 



AM 



1 A(r) 

c 2 r 



(B.8) 



(B.9) 



where A = Aj. This expression of the modification of -E suggests us to try a change of 
the electric dipolar moment p by which the modification produced by the substitution 
(IB.5j) can be compensated. Let be the substitution: 

p _> p' = p + Ap. (B.10) 
We try to determine Ap such that the modification ( |35i) of E is compensated. Since 



E(r,t) P ^ P E(r,t) 



E(r,t) 



1 



47T5 

1 

A'KSqC 2 



didj 



A Pj (r) 1 A^(r) 



V V 



r cr 
Ap(r)\ 



r 

1 Ap 

c 2 r 



(B.11) 
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the comparison with equation (1B.9P gives Ap = A/6c 2 and the substitution is in fact 



^ 6c 2 



P -> £> 

Let us consider also the change of £? by the substitutions ( IB. 51) and (IB. 121) 
B(r,t) 

= B{r,t) 
= B(r,t) 



(B.12) 



4-7T r 



/-to o Afc (r 



247TC 2 

/'■o 



Since 



24vrc 2 
* Afc(r 



r 

Ak (r) 



JM)_ 
24tt 

24tt 



e i^'fe djdid q 



" ~ 24^^ ^ ^ " 
(^ 9 A fc (T) + <^A 9 (r) + 5 gfe Aj(r) 



0. 



A fc (r) 
c 2 r 



and e ijk djd k = 0, 



£? is invariant to such a substitution. Concluding, we can substitute the tensor "P 1 ^ by 
the STF tensor if we perform simultaneously the substitution p — > p' = p + A/6c 2 . 

Since we have to consider, together with the contributions of the electric octupolar 
moment to the ^-singularities of the field, such contributions of the magnetic quadripolar 
moment also, we must show how one can substitute the momentum tensor M^ 2 ^ by the 
corresponding STF projection MS 2 \ Let be the components 

2 f 

M y (t) = - / d 3 xxi(r x J(r,t)) . 
6 Jv 

Concerning the magnetic quadrupole moment lw 2 % we have a simple procedure for 
obtaining the STF projection (up to a factor). Let us write the identity: 



M 



M 



where the first bracket represents the symmetric part of this tensor, and the second, the 
antisymmetric one. The symmetric part is, for this case (n = 2), a STF tensor Ad^. 
Therefore, 

1 

where 



M w = Mij + -e ijk N fc , 



Ni 



(B.13) 

e kij = ^- J [r x (r x J)] k d 3 x = — J [{r • J)x k - r 2 J k ] d 3 x. (B.14) 



The substitution M^ 2 ^ — > MS 2 ^ in equation (j55j) gives 



E (r,t) ™ B(r,f) 

= £(r,t) + 

= £(r,t) + 



1 PP P flfl N * (r) 

6j Sijk^lkq OjOl 



1 






167re c 2 


v r ) 


1 


v 









r 

7V(t) 

c 2 r 



(B.15) 



Singular behaviour of the electromagnetic field 



21 



where N = Njej. If we perform the substitution of the electric dipolar moment 
p — > p" = p + A'p, the electric field E is modified by an expression given by equation 
( IB.llf) (with A — > A'). It is easy to see that the compensation of the changes in the 
electric field produced by the substitutions M^ 2 ^ — > MS 2 ^ and p — > p" = p + A'p, is 
given by the choice: 



(B.16) 



Concerning the change of B by the substitutions p — > p + A'p, M^ 2 ^ — > jV^ 2 -*, we can 
write: 

P ^ + A> P M^- B{ t) /xo v x A^(r) 

v y 4tt r 



167T 



' 99 9 W _ aA N,(r) 



16nc 2 



16n 



V V- (V x 



r 

iV(r) 
r 



1 _ JV(t) 
-V x — — 



r 

= B(r,t), (B.17) 

i.e. the invariance of B. Finally, we can express the total effect of the two substitutions 
p(3) _> j>(3) anc j |\/|(2) _^ ^v|(2) represented by the transformation of the electric dipolar 
moment: 

1 



or 



where 



p(3)_p(3) 1 ■■ M ( 2 )^Ai< 2 > - 1 v 

oc z ocr 4cr 

p(3)^p(3) 1^,(3)^^(3) , 1 • 

p — ► — ► P = P~ —T, 

& 



T=-N-~A 

4 6 

is the electric dipolar toroidal moment [12 
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